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Abstract. Sexual size dimorphism (SSD), the difference in body size between males and females, is a widespread phe¬ 
nomenon and commonly attributed to variation in sex-specific patterns of selection. Using data from three populations of 
Polypedates megacephalus , we tested the hypotheses that sexual selection favoured large male body size and that fecundity 
selection favoured large female size. Females had a larger body size than males in all populations and thus exhibited a fe¬ 
male-biased sexual size dimorphism. Standardized selection gradients showed that sexual selection for large male size was 
weak at best and hence unlikely to be an important driving force for SSD. As fecundity selection favouring large females 
was weak as well, it could not explain the female-biased SSD patterns found. However, a faster growth rate in females fa¬ 
voured a larger size, providing an explanation for female-biased SSD. Our findings suggest that female-biased SSD does 
not result from a phenotypic selection of female by size in P. megacephalus. 
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Introduction 

Sexual size dimorphism (SSD) is a widespread phenome¬ 
non among animals (Darwin 1871). The occurrence of SSD 
has formerly been explained by several hypotheses. The 
sexual selection hypothesis predicts that intrasexual com¬ 
petition or intersexual mate choice favours a larger/smaller 
body size in one sex (Andersson 1994). Alternatively, SSD 
may arise through the reproductive biology of the females, 
whereby a strong fecundity selection favours females that 
have an increased reproductive output due to their being 
larger (Herczeg et al. 2010, Han & Hu 2013, Ma et al. 2015). 
However, apart from these two evolutionary hypotheses, it 
is possible that SSD might simply arise as a consequence of 
sexual differences in age, an age-specific growth rate, or bi¬ 
maturity between the sexes (for amphibians see Halliday 
& Verrell 1988, Kupfer 2007, Zhang & Lu 2013, Liao et al. 
2013, Liao et al. 2015, Liao et al. 2016) and ecology/niche di¬ 
vergence (Selander 1966, Serrano-Meneses et al. 2006). 

Thus far, most SSD studies have focused on interspecific 
patterns in amphibians (e.g.. Shine 1978,1979, Monnet & 
Cherry 2002, Liao et al. 2013, Zhao et al. 2016), where¬ 
as it was only during recent years that proximate mecha¬ 


nisms resulting in SSD at intraspecific level were analysed 
in wild populations (Liao & Chen 2012, Liao 2013, Liao 
et al. 2015). Consequently, most of these studies have fo¬ 
cused on taxa with female-biased SSD, which has indicated 
the positive selection on female body size to be stronger 
than that on male body size (Liao 2013, Liao et al. 2015). 
Females are larger than males in 90% of the species in an- 
urans (e.g.. Shine 1979, Kupfer 2007, Zhang & Lu 2013). 
In most species, a larger female produces a larger clutch 
size and thus has a higher fecundity. Female-biased SSD 
in anurans is thought to result from a larger clutch size/ 
fecundity, which favours a larger female size (Monnet & 
Cherry 2002, Han & Hu 2013). However, studies on phe¬ 
notypic selection on male and female sizes resulting in fe¬ 
male-biased SSD by comparing the standardized selection 
gradients are rare in amphibians (Liao et al. 2015). 

The aim of this study was to investigate proximate mech¬ 
anisms that may result in a female-biased SSD of the spot¬ 
legged treefrog (Polypedates megacephalus). We analysed 
male and female body sizes from three populations along 
an altitudinal gradient in Guizhou Province in China. We 
investigated whether female fecundity increases with in¬ 
creasing body size. We also compared the strength of selec- 
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tion acting on males and females for all populations with 
sufficient data on male mating success and female fecundi¬ 
ty. Finally, we identified the age of all individuals to control 
for age-related differences in body size and SSD using skel- 
etochronological methods. This allowed us to evaluate the 
observed SSD patterns as results of sex-specific differences 
in the age structure within populations. 

Materials and methods 

Study species 

The spot-legged treefrog (Polypedates megacephalus ) is a 
medium-sized anuran, which is widely distributed in Chi¬ 
na and lives at altitudes ranging from 520 to 2,200 m (Fei & 
Ye 2001). Breeding activities begin in mid-April and con¬ 
tinue until late July. It is a lekking species in which males 
congregate at ponds to wait for females, whereas females 
will only be present at these pools during individual nights 
for mating, and males will then be quick to approach ar¬ 
riving females (Chen et al. 2016). Amplectant pairs move 
about at the edge of the ponds and will then produce foam 
on the ponds’ water surface. This will then attract other 
males that come to join the pair and so form an amplecting 
group of two to five males and one female. In males, no di¬ 
rect physical competition for possession of the females can 
be observed during the spawning period. 

Study sites 

We visited three P. megacephalus populations located at dif¬ 
ferent altitudes in Guizhou Province in western China dur¬ 
ing the 2014 and 2015 breeding seasons. The low-altitude 
population was located at 449 m in the Fanjing moun¬ 
tains (io8°44.67’ E, 27°46.33’ N), where frogs reproduced 
in a farm pond (5.0x3.6x2.8 m; lxwxd) framed by a farm¬ 
house. The mid-altitude population was located at 680 m 
in Shangzhong (io8°43.38’ E, 27°23.38’ N), where the frogs 
likewise reproduced in a farm pond (3.0X4.2X3.8 m) near a 
small pig farm. The high-altitude population was located at 
1,300 m in the Leigong mountains (io8°io.27’ E, 26°22.73’ N), 
where frogs reproduced in two natural ponds (2.0X1.5X1.2 m 
and 12.0x4.0x0.3 m, respectively) framed by a farmhouse. 

Sampling 

Field observations were conducted at the three study sites 
during the breeding season in two successive years. Af¬ 
ter each female had deposited her eggs, we captured all 
males and females at the breeding sites and then estimat¬ 
ed the clutch size of each nest. We collected a total of 209 
tree frogs (Leigong: 28 unmated males, 20 mated males, 
and 20 females; Shangzhong: 27 unmated males, 21 mat¬ 
ed males, and 21 females; Fanjing: 20 unmated males, 21 
mated males, and 21 females) during the two breeding sea¬ 
sons. All individuals were confirmed to be adults by direct 


observation of secondary sexual traits (presence of nuptial 
thumb pads in adult males and the inflated abdomen car¬ 
rying developing eggs of adult females) and their breeding 
behaviour. Body size (snout-vent length: SVL in mm) of 
each individual was measured to the nearest 0.1 mm using 
dial callipers. The second phalange of the longest toe of the 
right hind limb of each treefrog was removed and stored in 
4% neutral buffered formalin for age identification (see be¬ 
low). All treefrogs were then individually released at their 
original capture sites. 

Age identification 

Age was identified by means of skeletochronology, which is 
based on counting the number of lines of arrested growth 
(LAGs) in stained cross sections of the phalangeal bones 
(Castanet & Smirina 1990). We used an improved meth¬ 
od of paraffin-sectioning and Harris’s haematoxylin-stain- 
ing to produce histological sections for aging the adult ani¬ 
mals. This method has been used in the skeletochronology 
of many Chinese anurans before (Liao & Lu 2010, Liao & 
Lu 2012, Li et al. 2013, Huang et al. 2013, Jin et al. 2016). We 
selected cross-sections (13 pm thick) of the phalanx with 
the smallest medullar cavity for mounting them on glass 
sides. Digitalised photographs of cross-sections were taken 
using a Motic BA300 digital camera mounted on a Moti- 
cam 2006 light microscope at 400X magnification. After 
this, the number of lines of arrested growth (LAG) was re¬ 
corded by the same researcher (Chen Cheng). We con¬ 
firmed LAG endosteal resorption and double lines. Of 209 
adult specimens, 207 (145 males and 62 females) exhibited 
clear LAGs in their bone cross sections. Growth was as¬ 
sessed by a non-linear regression in SPSS (21.0) statistical 
package (Statistical Product and Service Solutions Compa¬ 
ny, Chicago, USA) using von Bertalaneffy’s (1957) equa¬ 
tion of St = S (i-e~ kt+b ), where S is body size at age't', S 
is the estimated asymptotic maximum size, 'k' is a growth 
coefficient, and 'b' is a constant. The growth rate can then 
be calculated as R = dS/dt = k (S -S t ). 

v max v 

Clutch size 

A total of 62 amplectant treefrogs were collected in two 
breeding seasons and transported to the laboratories that 
were close to the breeding sites. The pairs were kept sepa¬ 
rately in tanks (20x20x30 cm) filled with pond water to 
allow oviposition. We collected a total of 54 clutches and 
counted the total number of eggs in each clutch in a dish. 
After the experiments, all individuals and eggs were re¬ 
turned to the places where they had been collected. 

Selection gradients 

To compare strengths of selective forces on male and fe¬ 
male size, linear standardized selection gradients were 
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Table 1. Mean body sizes (mm), age (years), and growth rates for three Polypedates megacephalus populations in Guizhou Province, 
China. (Specimens: Fanjin: males, n = 42; females, n = 21; Shangzhong: males, n = 48; females, n = 21; Leigong: males, n = 58; fe¬ 
males, n = 20). 


Population 

Body size 

Males 

Min. age Average age 

Max. age 

Growth rate 

Body size 

Females 
Min. age Average age 

Max. age 

Growth rate 

Fanjin 

43.6±3.7 

2 

2.6±2.0 

4 

1.47 

59.3±3.1 

2 

2.8±0.8 

4 

2.61 

Shangzhong 

43.0±2.9 

2 

2.9±0.8 

5 

1.31 

57.4±4.5 

2 

2.9±0.6 

4 

3.65 

Leigong 

44.0±3.2 

2 

3.1±0.8 

5 

1.94 

58.4±5.4 

2 

3.0±0.8 

4 

2.54 


estimated (Lande & Arnold 1983) for each of the three 
populations from the data on male mating status and fe¬ 
male fecundity. Following the methods of Lande & Ar¬ 
nold (1983), body size was standardized (within each sex) 
by subtracting the population’s mean from individual val¬ 
ues and dividing them by the samples standard deviation. 
Relative fitness was calculated by dividing the individual 
fitness (clutch size [females] or mating status [males]) by 
the mean fitness (i.e., mean clutch size or mating status) in 
the given population. We calculated the standardized lin¬ 
ear selection gradients using a regression of relative fitness 
on standardized body size. Univariate linear selection gra¬ 
dients (P) were estimated from the model, relative fitness 
(w’) = c + (standardized body size). We also calculated 
the non-linear selection gradients by applying a multiple 
regression of relative fitness on standardized body size and 
its square. Univariate non-linear selection gradients {y) 
were estimated from the model, relative fitness (w’) = c + P 
(standardized body size) + P 2 (standardized body size) 2 , 
where y = 2(U 

Statistical analyses 

All analyses were performed using the SPSS (21.0) statisti¬ 
cal package. Body size and clutch volume were log io -trans- 
formed. To test for sex differences in mean body size or 
age between the three populations, we used a Generalized 
Linear Model (GLM) with 'logj (SVL) or 'age' as the de¬ 
pendent variable, 'sex' as a fixed factor, and 'population' as a 
random effect. We used a Generalized Linear Mixed Mod¬ 
el (GLMM) with 'log io ' (SVL) as the dependent variable, 
'sex' as a fixed factor, and 'population' as a random effect, 
'age' and 'altitude' as covariates to test for sex differences 
in relative body size between populations. To assess the ef¬ 
fect of female body size on clutch volume, another GLMM 
was used. Here, 'population' was included in the model as 
a random effect, 'year of collection' as a fixed factor, and 
'female body size', 'female age', and 'altitude' as covariates. 
We tested for the variation of SSD among populations with 
'age' being added to the model as a covariate together with 
'sex x age' (fixed effect) and 'sex x population' (random ef¬ 
fect) interactions. A significant sex x age interaction would 
be suggestive of sex differences in growth rates, whereas a 
significant sex x population value would provide evidence 
for between-population variation in SSD. 


Results 

Females had a larger body size than males (GLM: F t = 
1009.008, P < 0.001; Table 1), but the average body size 
did not differ significantly between the three populations 
(F 2 9s = 2.604, P = 0.277). Likewise, there was a non-signifi¬ 
cant difference in average age between populations (F 2 2oj = 
7.280, P = 0.121) and between the sexes (F 2 2o3 = 0.204, P = 
0.696; Table 1). Females always had a larger body size than 
males (F | 203 006 = 70.602, P < 0.001) when considering the 
effect of age on body size (F 2 , D3l82 = 16.089, P < 0.001). Al¬ 
titude did not affect average body size between the three 
populations (GLMM: altitude: F 2 50409 = 2.108, P = 0.153; 
population: Z = 0.560, P = 0.576). The sex x age interaction 
revealed a non-significant difference in the relationship be¬ 
tween body size and age (=growth rate) between the sexes 
(F = 4.467, P = 0.018; Table 1). We also did not observe 
a significant sex x population variation between popula¬ 
tions (Z = 0.512, P = 0.609). 

The GLM revealed that clutch size differed significantly 
between populations (F 2 _ 3 = 95.859, P < 0.001), but not rel¬ 
ative to female body size (F 2 53 = 2.442, P = 0.097). Clutch 
size (286.4 ± 30.0) was highest in the high-altitude popula¬ 
tion and the natural ponds whereas it was lowest (192.4 ± 
21.8) in the farmed frogs at the lowest altitude (Fig. 1). 
Clutch size did not increase with increasing female body 
size (GLMM: Fj ;o oo = 1.039, P = 0.313; Fig. 1), even when 
the age effect on body size was considered (body size: 
F 0 = 1.171, P = 0.285; age: F , = 1.101, P = 0.299). 

1, 49.098 ' O 1, 49.000 

The standardized linear or non-linear selection gra¬ 
dients demonstrated that sexual selection did not favour 
large males in the three populations, and the overall selec¬ 
tion gradient estimated across all populations was non-sig¬ 
nificant. Meanwhile, the standardized linear or non-linear 
selection gradients on female size revealed that fecundity 
selection did not favour large females in all three popula¬ 
tions, and the overall selection gradient was non-signifi¬ 
cant as well (Table 2). 


Discussion 

Our results uncovered significant differences in body size 
between males and females (i.e., sexual size dimorphism: 
SSD) for all the spot-legged treefrog ( P. megacephalus) 
populations investigated, even after our correcting for age 
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effects. Furthermore, the differences in the growth rate 
caused by differences in altitude between populations and 
between the sexes exhibited a strong influence on body 
size. The patterns of SSD in the three populations can be 
explained as being a simple consequence of the popula- 
tional variation of the growth rate. The standardized se¬ 
lection gradients from wild populations did not provide 
strong evidence of a positive selection favouring large male 
size. Hence, sexual selection in favour of larger males can¬ 
not explain SSD in P. megacephalus. As fecundity selection 
did not favour larger females, it failed to explain the fe¬ 
male-biased SSD. This was supported by the discovery that 
female fecundity was not correlated with female body size 
both within and between populations. 

There is some evidence from interspecific studies of pat¬ 
terns of SSD in vertebrates (Fairbairn 1997), irrespective 
of whether females (as in most ectotherms) or males (as in 
most birds and mammals) are the larger sex (Szekely et al. 
2004), which needs to be looked at. The majority of studies 
suggest that sexual selection on male size is likely to explain 
male-biased SSD (Andersson 1994). In a few cases, fecun¬ 
dity selection has been suggested to explain female-biased 
SSD (Herczeg et al. 2010, Liao et al. 2015). Female-biased 
SSD is the prevalent pattern in anurans (Shine 1979), and 
this pattern was true also in the three P. megacephalus pop¬ 
ulations analysed. Selection gradients on both males and 
females could have produced SSD. In most cases, male- 
biased SSD is generally driven by the preferred selection 
on larger males as a consequence of sexual selection re¬ 
sulting from male-male competition or large-male-advan¬ 
tage (Fairbairn & Preziosi 1994). In this study, we found 
that the male selection gradients were non-significant for 
all populations, thus demonstrating that sexual selection 
did not favour larger males in P. megacephalus. In contrast, 
sexual selection might in fact have favoured smaller males 
in taxa with female-biased SSD, because of the increased 
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Figure 1. A non-significant relationship between female body size 
and clutch size in Polypedates megacephalus in the three populations 
investigated (female specimens: Fanjin, n = 19; Shangzhong, n = 17; 
Leigong, n = 18). Clutch size was largest in the Leigong Mountain 
population and smallest in the Fanjin Mountain population. 
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Table 2. Standardized selection gradients for sexual selection on 
male body size (dependent variable, mating status [mating or 
single]), and fecundity selection on female body size (depend¬ 
ent variable, reproductive fitness [clutch size/ mean clutch size]). 
Univariate linear selection gradients (|5) were estimated from 
the model, relative fitness (w’) = c + (5 (standardized body size); 
Univariate nonlinear selection gradients (y^ were estimated from 
the model, relative fitness (w’) = c + (3 (standardized body size) + 
(3, (standardized body size) 2 , where y = 2P 2 . 


Population 

Males 

Females 


Pi 

V, 

Pi 

Yi 

Fanjin 

-0.071 

0.022 

-0.003 

-0.004 

Shangzhong 

-0.002 

0.062 

-0.001 

0.002 

Leigong 

-0.079 

0.080 

-0.007 

-0.032 

Overall 

-0.052 

0.032 

-0.006 

0.006 


mobility or agility of small males (see Szekely et al. 2004 
for birds). However, we did not observe any markedly im¬ 
proved mobility or agility in male P. megacephalus. Alter¬ 
natively, fecundity selection has been suggested to explain 
female-biased SSD due to the increased reproductive in¬ 
vestment of large females (e.g., Herczeg et al. 2010, Liao 
et al. 2015). We found that the standardized estimates of 
the strength of selection on females did not exhibit a posi¬ 
tive selection, suggesting that fecundity selection did not 
favour large females, which is inconsistent with results that 
have previously been observed in other anurans ( Rhaco- 
phorus omeimontis: Ma et al. 2015; Bufo andrewsi : Liao et 
al. 2015). 

Intersexual differences in life-history traits, such as age, 
can also drive a female-biased SSD in anurans (Monnet & 
Cherry 2002, Zhang & Lu 2013). In P. megacephalus, the 
higher average female age and longevity did not contribute 
to a female-biased SSD. However, previous studies in an¬ 
urans suggested that age and growth rates can significantly 
influence body sizes within each sex (Shine 1979, Rein- 
hard et al. 2015). Consequently, differences in the growth 
rate and/or the time available for growth, as well as the two 
factors acting simultaneously, can result in SSD. As females 
have experienced a longer duration of growth than males, 
SSD will increasingly be biased towards females (Zhang 
& Lu 2013). In this study, we found that males and females 
did not exhibit differences in age at sexual maturity and 
longevity, but there were significant differences in growth 
rate (as an indicator of the sex x age interaction on body 
size), demonstrating that females attained larger body siz¬ 
es and grew faster towards the asymptotic body size than 
their male conspecifics. Consequently, faster growth re¬ 
sults in more energy being devoted to somatic growth in 
order to achieve a larger body size, which leads to a female- 
biased SSD in this species, namely all amphibians faced the 
trade-off between reproduction and growth. Moreover, 
male and female growth are related to age and sexual ma¬ 
turity. As a result, males may even have the same growth 
rate, but may mature earlier than females, thus resulting in 
a smaller male size in some species (Duellman & Trueb 
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1986, Morrison & Hero 2003). However, we found sex 
x population to be non-significant between populations, 
suggesting that phenotypic selection on certain females did 
not result in variation in female-biased SSD selection. 

In summary, our study suggests that sexual selection 
does not favour male body size in P. megacephalus, and 
cannot explain the occurrence of the female-biased SSD 
either. Meanwhile, fecundity selection does not favour a 
large female size, which is also unlikely to explain the fe¬ 
male-biased SSD and an inference supported by the find¬ 
ings that female fecundity is not significantly correlated 
with female body size in P. megacephalus. However, faster 
female growth rates indicate that females apply more en¬ 
ergy to growth and grow larger as a result (Morrison & 
Hero 2003). Moreover, the predation risk in terms of di¬ 
versity and density of predators in the different habitats 
(i.e., farm and natural ponds) at different altitudes impacts 
on growth rates in anurans (Beck et al. 2008). Hence, dif¬ 
ferences in growth rates (as an indicator of the sex x age 
interaction on body size) between the sexes is a likely ex¬ 
planation for the female-biased SSD. 
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